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Abstract The phase-space structure of two families of galactic potentials is 
\—^ | approximated with a resonant detuned normal form. The normal form series 

is obtained by a Lie transform of the series expansion around the minimum of 
the original Hamiltonian. Attention is focused on the quantitative predictive 
ability of the normal form. We find analytical expressions for bifurcations 
of periodic orbits and compare them with other analytical approaches and 
with numerical results. The predictions are quite reliable even outside the 
convergence radius of the perturbation and we analyze this result using re- 
summation techniques of asymptotic series. 

^ ! 
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1 Introduction 

o 



Normal forms are of invaluable help in approximating the phase-space struc- 
ture of non-integrable systems and in providing a picture of the dynamics 
in the regular domains (Giorgilli and Locatelli, 2006). In particular, detuned 
resonant normal forms allow us to investigate in detail several features of 
non-linear oscillators. Leaving free parameters, these tools have been exten- 
sively applied to get useful qualitative information (Cushman and Bates, 
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1997; Broer at al., 2003). However, quantitative predictions in concrete ap- 
plications are as much as useful. To take a specific example, let us consider 
the bifurcations of periodic orbits: the stability-instability threshold in a 
suitable parameter space can be determined by studying the nature of criti- 
cal points of cither exact invariant functions in the normalizing variables or 
of approximate integrals in the original variables. The agreement between 
both methods is obtained through series expansions in characteristic param- 
eters (Belmonte et al., 2007), for example the energy and the ellipticity of 
the equipotentials: in the case of the galactic 'logarithmic' potential (Binney 
and Tremaine, 1987), the agreement with other analytical or numerical ap- 
proaches is very good even truncating the normal form at a relatively low or- 
der. Attempts have also been made to treat perturbations of the isochronous 
sphere (Gerhard and Saha, 1991; Yanguas, 2001). 

Aim of this work is to study in a more systematic way the quantitative 
predictive ability of these expansions. We investigate the result of different 
choices of the effective perturbation order of the detuning term. The usual 
approach to construct detuned normal forms is based on the idea of consider- 
ing detuning zero-order terms as being of higher order and therefore treating 
them as part of the perturbation. In the process of normalization, new terms 
appear at orders (and with coefficients) that depend on the 'effective' order 
of the detuning term. In turn, different terms appear in the process of inver- 
sion to construct approximate 'third' integrals of motion (Contopoulos, 2004; 
Giorgilli, 2002). Benchmarks used to test the predictions are in all cases the 
outcome of other analytical approaches. The results point out that the best 
choice is always that of treating the detuning term as the lowest non-zero 
order term compatible with the symmetries of the problem at hand. 

A remarkable side-effect of expressing the stability-instability threshold 
as a series expansion, is that its predictive ability goes well beyond the radius 
of convergence of the perturbing expansion. In the 'galactic' cases studied 
here, we show that the bifurcation energy of the major-axis periodic orbit is 
predicted with high accuracy up to values much higher than the 'harmonic 
core' energy. We are using probably divergent series; however, their trunca- 
tions stay close to actual functions. To justify these extended results, we try 
to estimate an optimal truncation by exploiting some resummation technique 
based on asymptotic series and continued fractions. 

The plan of the paper is as follows: in Section 2 we recall the procedure 
of detuned resonant normalization as applied to double-reflection symmetry 
potentials; in Section 3 we sketch the basics of the theory of stability of 
normal modes for these systems; in Section 4 we study the effective order 
of the detuning and compare the predictions with numerical data, analyzing 
the results using resummation techniques of asymptotic series; in Section 5 
we give our conclusions. 

2 The Hamiltonian system 

We are interested in 2-degrees of freedom natural systems of the form 

1 
2 



H(p,d) = -{pl+p 2 y ) + V(x 2 ,y 2 ), (1) 



with V a smooth potential with an absolute minimum and reflection symme- 
try with respect to both axes. The motivation for the choice of this symmetry 
stems from its interest in problems of galactic dynamics. In particular, we 
will examine the following two examples 

V L = Uog{R 2 + x 2 + y 2 /q 2 ), (2) 



V C = y/R 2 + x 2 + y 2 /q 2 - R. (3) 

Both choices are useful models for elliptical galaxies: in particular, the loga- 
rithmic potential Vl is used to model a core embedded in a dark matter halo. 
For every finite values of the "core radius" R, the choice R = 1 can be done 
without any loss of generality. However, in both cases it is also of relevance 
the singular limit R — > associated to central density cusps. At high energy, 
the dynamics in the non-singular cases tend to the corresponding scale-free 
singular limits, respectively the singular logarithmic potential Vls and the 
"conical" potential Vcs ■ Since it is difficult to devise a perturbative approach 
to approximate the dynamics in the singular cases, one may hope to infer 
information on them from the study of the non-singular cases at sufficiently 
high values of the energy. 

With the choice R = 1, the energy E may take in both cases any non- 
negative value 

< E < oo. (4) 

The parameter q gives the "ellipticity" of the figure and ranges in the interval 

0.6<q<1. (5) 

Lower values of q can in principle be considered but correspond to an un- 
physical density distribution. Values greater than unity are included in the 
treatment by reversing the role of the coordinate axes. In order to implement 
the normalization algorithm, the Hamiltonian has to be expressed as a series 
expansion around the equilibrium. Performing the scaling transformation 

Py — > VQPy, V — ► —^ ( 6 ) 

the Hamiltonian can be expanded as 

where 



q 

k=0 y ' k=0 



s = x 2 + ly 2 (8) 

and bk are real coefficients (with bo — 1/2). The series expansions of the 
potentials (|2I3[) around the origin, with R= 1, are: 

2 2-22-3 2-4 v; 



and 

1 1-1 o 1-1-3 , 1-1-3-5 , 

Vc = r-2-4 S+ 2^-6 S -2^6~8 S+ --- (10) 

respectively. 

We look for a new Hamiltonian in the new canonical variables P, Q, given 

by 

oo 

tf(P,Q) = ^if„(P,Q), (11) 

with the prescription that 

{H ,K} = 0. (12) 

The zero order (unperturbed) Hamiltonian, 

H Q = ^(P^+X 2 ) + j-(P 2 +Y 2 ), (13) 

with unperturbed frequencies u>\ = 1 and io 2 = 1/q, is expressed in terms of 
the new variables found at each step of the normalizing transformation. In 
these and subsequent formulas we adopt the convention of labeling the first 
term in the expansion with the index zero: in general, the 'zero order' terms 
are quadratic homogeneous polynomials and terms of order n arc polynomials 
of degree n + 2. 

It is customary to refer to the normal form constructed in this case as a 
"Birkhoff" normal form (Birkhoff, 1927). The presence of terms with small 
denominators in the expansion, forbids in general its convergence. It is there- 
fore more effective to work since the start with a resonant normal form 
(Sanders & Verhulst, 1985), which is still non-convergent, but has the advan- 
tage of avoiding the small divisors associated to a particular resonance. To 
catch the main features of the orbital structure, we therefore approximate 
the frequencies with a rational number plus a small "detuning" 

— = q = — - + 5. (14) 

uo 2 m 2 

We speak of a detuned (nii:m 2 ) resonance, with mi + rri2 the order of the 
resonance and, performing the rescaling 

m 2 qH, (15) 
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redefine the Hamiltonian in the form 
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-y 2 )}+m 2 {±5(pl+x 2 ) + qY,b k s k+1 }. 

k=Q " fc=l 

(16) 
The procedure is now that of an ordinary resonant "Birkhoff-Gustavson" 
normalization (Gustavson, 1966; Moser, 1968) with two variants: the coor- 
dinate transformations are performed through the Lie transform and the 
detuning quadratic term is treated as a term of higher order and put in the 
perturbation. 



The new coordinates P, Q result from the canonical transformation 

(p,q)=T x (P,Q). (17) 

Considering a generating function x, the Lie transform operator T x is defined 
by (Boccaletti and Pucacco, 1999) 

oo 

T X ^J2 M >< ( 18 ) 

fc=0 

where 

k . 

M = 1 , M k = J2 rL Xj M k -j . (19) 

3=1 

The functions Xfe are the coefficients in the expansion of the generating func- 
tion of the canonical transformation and the linear differential operator L g 
is defined through the Poisson bracket, L g (-) = {g, ■}. 

The terms in the hew Hamiltonian arc determined through the recursive 
set of linear partial differential equations (Giorgilli, 2002) 

K = H Q , 

K\ = TCi + MiHq = Tii + L Xl Tio, 



Ki = H.2 + M{H\ + M2TC0 = TL2 + L Xl TLi + 2^x1^° "*~ ^x^^-o, 



(20) 



K n = H n + E"=i M^jHj + M n H , 

where L"z(-) — {<?, {g, •}}. 'Solving' the equation at the n-th. step consists of a 
twofold task: to find K n and Xn- T ne unperturbed part of the Hamiltonian, 
Ho, determines the specific form of the transformation. In fact, the new 
Hamiltonian K is said to be in normal form if, analogously to JT3 



{H o ,K} = 0, (21) 

is satisfied. We observe that, in view of the reflection symmetries of the 
potentials (|2I3[) . in the chain (|2H)) they appear only terms with even index 
and so the normal form itself is composed by even index terms only. 

We have to discuss how to treat the detuning term: it is considered as a 
higher order term and the most natural choice is to put it into H.2- However, 
there is no strict rule for this and one may ask which is the most 'useful' 
choice, always considering that applications are based on series expansions 
with coefficients depending on q. We remark that, different choices of the 
effective order, say d, of the detuning, lead to different terms of higher order 
in the normal form. We also observe that, whatever the choice made, the 
algorithm devised to treat, step by step, the system (|2"U)) must be suitably 
adapted to manage with polynomials of several different orders. In practice, 



since at each step the actual order of terms associated to detuning is lower 
than the corresponding effective order, the algorithm is adapted by incorpo- 
rating routines already used at previous steps. In practice, at step say j, we 
have an equation of the form 

K t = Hi + A, + 8B^ d + 6 2 B 3 _ 2d + ... + L Xj H , (22) 

where Aj, Bi are homogeneous polynomials of degree i + 2 coming from pre- 
vious steps. As usual, the algorithm is designed to identify in all terms with 
the exclusion of 

L x .H = -L no ( Xj ), (23) 

monomials in the kernel of the linear operator L-h . These monomials are used 
to construct Kf the remaining terms are used to find Xj hi the standard way. 
It is clear that both the normal form and the generating function are affected 
by the effective order of the detuning term. 

In both cases (|9ll0p investigated here, with the detuning treated as a 
term of order 2, the next appearance of a related term is in Kq. Rather, if it 
is treated as a term of order 4, the next appearance of a related term is in 
K%. Truncating at order 6 (polynomials of degree 8) is therefore sufficient to 
make a comparison with other predictions not sensitive to the detuning. 



3 Stability of periodic orbits 

A very useful setting in which to perform quantitative predictions is that 
of the stability threshold of normal modes and/or periodic orbits in general 
position (de Zeeuw and Merritt, 1983; Fridman and Merritt, 1997). In Bel- 
monte et al. (2007) we have shown that these predictions can be obtained 
both exploiting the nature of critical points of the normal form constrained 
on the manifold determined by the first integral and studying fixed points on 
a surface of section constructed using the approximate integral in the original 
coordinates. The agreement between the two methods is up to the level of 
the truncation if the transition curve in the parameter space is expressed as 
a series expansion around the given resonance. The choice of the method can 
be done simply on the basis of computational simplicity: for example, if one is 
interested in the instability transition of the normal modes, to work with the 
normal form is usually easier. Moreover, it is simpler to use 'action-angle-Zifce' 
variables, defined through the transformation 

X = v^ cos 0i, P x = i/2Jisin0i, (24) 

Y = v/2J2Cos0 2 , Py = v / 2^2sin0 2 . (25) 

A case that is both representative of the state of affairs and useful in galactic 
applications is that of the stability of the a;-axis periodic orbit (the 'major- 
axis orbit', if q is in the range |5j)). Among possible bifurcations from it, 
the most prominent is usually that due to the 1:2 resonance, producing the 
'banana' and 'anti-banana' orbits (Miralda-Escude and Schwarzschild, 1989). 
We will mostly investigate this problem in the potentials (12|3j) and will briefly 
discuss other less relevant cases. 
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Let us consider the l:2-symmetric resonant normal form, so that the 
frequency ratio (|14p now is 

uj 1 /uj 2 = 1/2 + 5. (26) 

The lowest order incorporating the resonance is 4 and we can write the 
corresponding terms as 

K a =H = Ji + 2J 2 , (27) 

K 2 = 25J 1 -P^{J 1 ,J 2 ), (28) 

K 4 = P (3) ( Ji, J 2 ) + kjf J 2 008(40! - 20 2 ). (29) 

where the polynomials P^ and P^ are homogeneous of degree 2 and 3. In 
the logarithmic and conical case they are respectively 

and 

p|» - 1 ^ + -H) + 2 J ' J - 

The constant fc in front of the resonant term in K4 is respectively q/8 in 
the logarithmic and (1 + (?)/32 in the conical case. For simplicity we have 
included the detuning term, 2(5 Ji, in K 2 . The generating functions have, in 
both cases, xo — 1 and xi = 0. The first non trivial term is 

X2 = -- J1M2 + Ji) sin26>i - ■£- J\ sin46>i + 
2 lb 

-7J2 I Ji + —J sin 20 2 - -Lj2 s in40 2 + 

+ ^ Ji -h ( sin(26»i - 20 2 ) - \ sin(20! + 20 2 ) 

in the logarithmic case and 2^2 in the conical case. 

Action-angle-like variables are singular on axial orbits. We have to use 
mixed variables: action-angle variables on the normal mode and Cartesian 
variables on the normal bundle to it. To analyze the stability, we determine 
the condition for the normal mode to be a critical curve of the Hamiltonian 
in these coordinates and assess its nature by considering the function 

K M =K + f iH 0l (30) 

where /1 has to be considered as a Lagrange multiplier to take into account 
that there is the constraint 

Wo = e (31) 



associated to the existence of the second integral. The Lagrange multiplier is 
found by imposing 

dif (/i) =0, (32) 

that is the total differential of ((30}) vanishes on the normal mode. Its nature 
is assessed by computing the matrix of second derivatives of K^>: if the 
Hessian determinant of the second variation is positive definite, the mode 
is elliptic stable; if it is negative definite, the mode is hyperbolic unstable 
(Kummer, 1977; Contopoulos, 1978). 

The equation det[d 2 K^(£)] — is an algebraic equation of degree M in 
£, where M is the order of truncation of the normal form. At order M = 4 
it is 

(48 - 96q + 24(3<? - l)£ + (26 - 153g + 153q 2 )£ 2 ) 

x (48 - 96q + 24(3g - l)£ + (26 - 159<? + 153q 2 )£ 2 ) = (33) 

in the logarithmic case (Belmonte et al. 2007) and, analogously, it can be 
shown to be 

(192 - 384 9 + 48(3g - l)£ + (41 - 213g + lbiq 2 )£ 2 ) 

x (192 - 384g + 48(3<? - l)£ + (41 - 219q + 153<? 2 )£ 2 ) = (34) 

in the conical case. 

The solutions of these equations give a set of curves in the q — £ plane. 
These are the loci of transition to instability of the normal mode and of 
bifurcation of a new family: among the set, the meaningful ones in the present 
case are those corresponding to the detuned 1:2 resonance. Therefore, the 
relevant solutions are those possessing a leading order going as 

£crit(<?) ~ Q~h (35) 

and this suggests to expand the selected solution as a power series in the 
detuning. However, in order to get a form usable in comparison with other 
results (for example coming from a numerical treatment) it is necessary to use 
a 'physical' energy variable rather than the parameter £. The conversion is 
possible if the physical energy E appears explicitly. According to the rescaling 
(|15p . we assume that m2qE is the constant 'energy' value assumed by the 
truncated Hamiltonian K. In the present instance rri2 = 2 so that, on the 
x-axis orbit, the new Hamiltonian is a series of the form 

K = 2q£-jq£ 2 + ... = 2qE. (36) 

The series (|3"6l can be inverted to give 



£ = E+^E 2 + ... (37) 



and this can be used in the treatment of stability to replace £ with E. The 
solutions can therefore be expressed as 

M/2 k 

fc=i ^ ' 





Potential Vl 


Potential Vc 
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Banana 


Anti — banana 


Banana 


Anti — banana 


1 


8 


8 


16 


16 


2 


20 
3 


28 
3 


248 
3 


536 
3 


3 


268 
9 


460 
9 


3608 
9 


18584 
9 


4 


1724 

27 


3928 

27 


43328 

27 


657848 

27 


5 


79184 


267404 


525704 


23668304 


405 


405 


81 


81 


6 


567178 
1215 


510200857 
405 


28118794 
1215 


4304374384 
1215 


7 


30991946 
25515 


615376795556 
8505 


309430864 
3645 


31575390356 
729 



Table 1 Coefficients in the expansion (|38p with M = 14 for the logarithmic poten- 
tial (banana, 2nd column and anti-banana, 3rd column) and the conical potential 
(banana, 4th column and anti-banana, 5th column). 



and in this form they can be used for quantitative predictions. The reason 
for a truncation order M/2 for this series is due to the fact that we now 
work with the energy rather than with phase-space coordinates: since M is 
even, the order is always an integer. A similar procedure can be followed 
for other bifurcations, keeping in mind that the lowest order to be included 
in the normal form in order to capture an m\:m<i symmetric resonance is 
2 x (mi + m 2 — 1) (Tuwankotta and Verhulst, 2000). Below, we will briefly 
touch on the case of the 1:1 resonance leading to the bifurcation of the loop 
orbit. 



4 Results and comparison with other methods 

In Belmonte et al. (2007) the series (|38f for the stability-instability transition 
of the resonant periodic orbits from normal modes have been computed for 
the potential ([9]) up to order M = 6. Here we extend those results to higher 
orders and to the other case of potential (jTTJ|> . 



4.1 Choice of the effective order of the detuning 



In Table [TJ we list the coefficients of (f38|) giving the bifurcation of the 1:2 
resonant periodic orbits ('banana' and 'anti-banana') for the logarithmic po- 
tential (JHJ) and the conical potential (fTTj|) . They have been obtained with a 
normal form truncated at order M = 14 and with the detuning treated as a 
term of order 2. There is a complete agreement with the analytical approach 
based on the Poincare-Lindstedt method (Scuflaire, 1995) and, as discussed 
below, there is a striking agreement with the numerical approach based on 
the Floquet method. 
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The results obtained by Scuflaire (1995) are based on treating the tran- 
sition as a parametric instability phenomenon and are therefore rooted in a 
quite different theoretical framework. He gets series relating the critical q to 
the amplitude on the axial orbit: once translated to the form (J3"B")) , the agree- 
ment of all fractional coefficients is complete up to the order we have arrived 
at, M/2 — 7. On the other hand, if the detuning is treated as a term of order 
d = 4 or greater, we get a disagreement in the coefficients starting from C3. 
This result confirms the analysis made above on the 'propagation' of the de- 
tuned terms in the normal form and show that the choice d — 2 is the optimal 
one. This choice has therefore been adopted for all subsequent predictions, 
which are compared with the numerical ones with a better agreement as long 
as the truncation order is increased. 

What is remarkable in all these results is that the numerical analysis are 
performed with the exact logarithmic (or conical) potentials (|2|3j) . whereas 
the analytical predictions are, in any case, produced with methods based 
on the series expansions (|9I10[) of the potentials with limited convergence 
radii. The usual attitude in normal form theory is to use the formal series as 
asymptotic series (Verhulst, 1996; Contopoulos et al., 2003; Efthymiopoulos 
et al., 2004) useful to approximate the dynamics in regular regions of phase 
space. We have adopted the same attitude in evaluating the reliability of the 
series giving the instability thresholds. 



4.2 Asymptotic series and resummation techniques 

Suppose we want to approximate a function f(x) with an infinite power series 
and define the "error" by truncating the series at order N: 

N 

e N (\x - sol) = f(x) - Y, c n{x - x ) n . (39) 

n=0 

We say that the series is asymptotic to the function if 

£n{\% — %o\) <C (x — x ) N , x — > xq, N fixed. (40) 

Typically the terms in the series get smaller for a while, but eventually they 
start to increase (the series diverge!). Since cn+i(x — xq) n+1 is an estimate of 
the error, we can find the optimal order of truncation N = N opt determining 
the smallest term (Bender & Orszag, 1978). The optimal order depends on 
the interval \x— Xo\' the larger the interval, the smaller iV op t and the accuracy 
in the approximation. On the contrary, for a convergent series, for every x 
within the interval \x — Xo\ 

00 
eiv(|a;-a;o|)= ^ c„(x-x )"^0, N -> 00. (41) 

n=N+l 

Once reached the optimal order, it can be disappointing to discard terms 
coming from a costly high-order computation. There are however other so- 
phisticated rules for 'summing' divergent series which make use of all terms 
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g 


N 


0.6 


0.7 


0.8 


0.9 


1 


0.800000 


1.60000 


2.40000 


3.20000 


2 


0.733333 


1.33333 


1.80000 


2.13333 


3 


0.763111 


1.57156 


2.60400 


4.03911 


4 


0.756726 


1.46939 


2.08680 





5 


0.758681 


1.53196 


2.56190 





6 


0.758214 


1.50208 


2.22160 





7 


0.758336 


1.51763 


2.48724 





E B 


0.758 


1.513 


2.401 


3.646 



Table 2 Subsequent truncations of expansion (|38[) with 
potential (banana). Eb is the value obtained by means 



M = 14 for the logarithmic 
of the Floquet method. 



(Bender & Orszag, 1978), like the construction of Pade approximant. A re- 
lated approach is that of constructing continued fractions from the original 
power series in the form 



m 



«o 



i + 



ai(x-xo) 



(42) 



i+- 



<■()> 



The coefficients a n can be computed from the c„ of the original series by 
expanding the continued fraction, comparing it with the original series and 
equating terms of the same order. Successive approximants obtained by trun- 
cating the fraction at various order may give an improvement in the asymp- 
totic convergence with respect to the original series (Khovanskii, 1963). In 
the following subsections we exploit these techniques to try to evaluate the 
optimal order of the series giving the instability thresholds. 



4.3 Bifurcation in the logarithmic potential 

We start by investigating the bifurcation of the 1:2 resonant periodic orbits 
('banana' and 'anti-banana') from the a;-axis orbit of the logarithmic poten- 
tial. Truncating to order 14, it is possible to get the loci in the q — E plane 
as in (|38|) with the coefficients of Table [TJ To evaluate these predictions, 
we treat the series (|38j) as an asymptotic series and evaluate truncations by 
computing the successive partial sums 



JY 



fc=l V Z 



N = l,...,M/2. 



(43) 



In Table [21 we report these partial sums for the banana, with q in the range 
from 0.6 to 0.9 and compare them with the numerical values obtained by 
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g 


N 


0.6 


0.7 


0.8 


0.9 


1 


0.800000 


1.60000 


2.40000 


3.20000 


2 


0.738462 


1.37143 


1.92000 


2.40000 


3 


0.753917 


1.45915 


2.14359 


2.81722 


4 


0.758120 


1.50678 


2.33218 


3.31985 


5 


0.758295 


1.51190 


2.37192 


3.51293 


6 


0.758323 


1.51399 


2.40779 





7 


0.758313 


1.51281 


2.38179 





E B 


0.758 


1.513 


2.401 


3.646 



Table 3 Subsequent truncations of the continued fraction (|42|) with M = 14 for 
the logarithmic potential (banana). 



means of the Floquet method (Miralda-Escude and Schwarzschild, 1989; Bel- 
monte et al. 2007) given in the last row. The numerical values of the partial 
sums are given with 6 digits just to show more clearly the asymptotic be- 
haviour: we can see that, up to q = 0.8, the predictions are apparently still 
(slowly) converging at N — 7. Only at the rather extreme value q = 0.9 we 
get an 'optimal' truncation order iVopt = 3, with a 10% error on the exact 
value of the critical energy. 

We may wonder if the continued fraction may help in speeding up the 
convergence rate: that actually this is the case can be seen in Table [3] where 
we report the partial sums computed with a continued fraction (|42p truncated 
at N G [1,7]. For all values of q up to 0.8, N ~ 6 is enough to reach a precision 
comparable to the numerical error. For q = 0.9 we get an optimal truncation 
order N opt = 5, with a 3% error on the exact value of the critical energy. 

Another check is that based on the more uncertain anti-banana transition, 
for which reliable numerical predictions exist only in the range 0.6 -r- 0.7. We 
put the results together in the same Table[U Even here, the continued fraction 
is clearly more efficient. 

Concerning the bifurcation of the loop orbit from the y-axis orbit, this 
happens in general at lower energies than that of the banana, therefore we 
presume to be predicted even better. This can be verified in the data of Table 
[5) where the continued fraction has been used. All partial sums are displayed 
since in any cases the optimal order is still unreached at N = 7: when entries 
are repeated is because the rounding off hides a truncation error which is 
still decreasing. 



4.4 Bifurcation in the conical potential 

In Table [6j we report the partial sums for the banana transition in potential 
([3]), with q in the range from 0.6 to 0.7 and compare them with the numerical 
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N 


0.6 


0.6 


0.7 


0.7 


1 


0.800000 


0.800000 


1.60000 


1.60000 


2 


0.893333 


0.905660 


1.97333 


2.08696 


3 


0.944444 


1.006320 


2.38222 


-2.32000 


4 


0.958993 


0.966571 


2.61499 


3.09306 


5 


0.965595 


0.969092 


2.82627 


3.38465 


6 


0.968026 


0.969675 


2.98186 


3.59815 


7 


0.969089 


0.969877 


3.11793 


4.15(*) 


AU 


0.970 


0.970 


4.292 


4.292 



Table 4 Subseque nt t runcations of series ((38} (first and third columns) and of 
continued fraction Q42p (second and fourth columns) with M — 14 for the loga- 
rithmic potential (anti-banana: Ea is the value obtained by means of the Floquet 
method.). The asterisk denotes a partial sum beyond the optimal truncation. 





q 


N 


0.6 


0.7 


0.8 


0.9 


1 
2 
3 
4 
5 
6 
7 


0.80000 
1.00000 
1.04000 
1.03738 
1.03776 
1.03783 
1.03783 


0.600000 
0.705882 
0.720000 
0.719333 
0.719408 
0.719419 
0.719419 


0.400000 
0.444444 
0.448000 
0.447891 
0.447900 
0.447901 
0.447901 


0.200000 
0.210526 
0.210909 
0.210903 
0.210904 
0.210904 
0.210904 


El 


1.038 


0.719 


0.448 


0.211 



Table 5 Subsequent truncations of continued fraction (|42[1 with M = 14 for the 
logarithmic potential (loop). El is the critical value obtained by means of the 
Floquet method. 



values obtained by means of the Floquet method. In Table [7J are reported 
the data for the loop orbit: even here, all partial sums are displayed since 
in any cases the optimal order is still unreached at N = 7: when entries are 
repeated is because the rounding off hides a truncation error which is still 
decreasing. 

In this potential, at higher values of the ellipticity, the bifurcation of 
the banana occurs, if any, at extremely high values of the energy: the same 
happens with the bifurcation of the loop at any q, with a critical energy much 
higher than in the logarithmic case. In both tables, it is clearly appreciable 
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g 


N 


0.6 


0.6 


0.65 


0.65 


0.7 


0.7 


1 


1.60000 


1.60000 


2.40000 


2.40000 


3.20000 


3.20000 


2 


2.42667 


3.31034 


4.26000 


10.6667 


6.50667 


-96.0000 


3 


2.82756 


3.20501 


5.61300 


9.22367 


9.71378 


113.029 


4 


2.98803 


3.80153 


6.42540 


13.0622 


12.2814 


-41.9262 


5 


3.05293 


3.09696 


6.91825 


7.67725 


14.3582 


23.1215 


6 


3.07607 


3.08855 


7.18186 


7.45669 


15.8394 


18.5519 


7 


3.08456 


3.08865 


7.32691 


7.46075 


16.9260 


18.6518 


Mb 


3.09 


3.09 


7.46 


7.46 


18.60 


18.60 



Table 6 Subsequent tru ncations of series (|38|l (first, third and fifth columns) and 
of continued fraction (|42|) (second, fourth and sixth columns) with M — 14 for the 
conical potential (banana). 





Q 


N 


0.6 


0.7 


0.8 


0.9 


1 


1.60000 


1.20000 


0.80000 


0.400000 


2 


oo 


4.80000 


1.60000 


0.533333 


3 


21.6000 


4.10323 


1.54074 


0.529870 


4 


8.50000 


3.64608 


1.51495 


0.529213 


5 


10.0922 


3.71312 


1.51728 


0.529239 


6 


9.92646 


3.70788 


1.51715 


0.529238 


7 


9.97972 


3.70929 


1.51717 


0.529238 


E t 


10.38 


3.710 


1.517 


0.529 



Table 7 Subsequent truncations of continued fraction (|42[) with M — 14 for the 
conical potential (loop). 



the accelerated convergence of the continued fraction up to those extreme 
values of the transition energy. 



4.5 Predictive ability of detuned normal forms 



The convergence radius of both expansions (j MTO)) in terms of the variable 
s is of the order unity. The corresponding energy level is E ~ 0.35 for the 
logarithmic potential Vl and E ~ 0.4 for the conical potential Vc- From 
the results presented above, it seems that the predictions obtained from the 
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normal form are quite reliable up to energy levels much higher than these. 
The two potentials examined here show a common orbit structure: however, 
features corresponding to the same phenomenon happen to occur at higher 
energies for Vc than for Vl- For example, from Tables [3] and [5] we see that, 
for q = 0.6 (namely, the ellipticity with the lowest values of the bifurcation 
energy), the banana bifurcates at E = 0.758 in Vl and E = 3.09 in Vc, 
with an astounding convergence of the analytical predictions (in particular, 
with the continued fraction). We are led to speculate about the possibility 
of extending this predictive ability to general features of the systems. More- 
over, we remark that at the energy levels reported here, both Vl and Vc 
slightly depart from the behavior of their scale- free counterparts. In view 
of the relevance of these results, it is of great interest the evaluation of an 
effective region of validity of the normal form. Therefore, we plan to extend 
the analysis to the study of the conservations of the approximate integrals of 
motion and to the explicit solutions of the equations of motion. 

5 Conclusions 

The results presented in this work provide a comprehensive setting to under- 
stand some stimulating but unsettled body of results obtained in previous 
studies. 

In Belmonte et al. (2006) we have started to investigate the stability of 
axial orbits in the logarithmic potential, using a normal form truncated to the 
first order incorporating the resonance: we got quite a satisfactory agreement 
with other analytical and numerical predictions but, among other things, we 
pointed out the troubles due to sensible differences between results coming 
either from the normal form itself ('final' normalizing variables) or from the 
approximate integral ('initial' physical variables). 

In Belmonte et al. (2007) we made a substantial step forward based on 
higher order normal forms: we were able to reconcile the results in the two 
different sets of variables by presenting the predictions as suitable power 
series in the same parameter space and, truncating at least at order M = 6, 
we provided good quantitative predictions. 

In the present paper we have tried to answer the natural question about 
the limits of validity of those predictions. Adopting the same approach usu- 
ally followed in exploiting approximate invariants constructed in a high order 
perturbation theory, we have shown how the power series representing the 
instability thresholds in the relevant parameter space can be interpreted as 
asymptotic series. As such, their truncated sums can be used to find an opti- 
mal truncation order and, in case, they can be resummed with the technique 
of the continued fraction, obtaining an improvement in the convergence rate. 
The generality of this setting allows us to conjecture that those results can be 
extended to arbitrary resonances and to periodic orbits in general position. 
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